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SUMMARY 


The  original  goal  of  the  present  contract  was  (i)  to  study  the  properties  of  mode 
conversion  from  magnetostatic  forward  volume  waves  (MSFVW's)  to  magnetostatic 
backward  volume  waves  (MSBVW’s)  through  a  region  of  bias-field  discontinuity,  and 
then  (ii)  to  investigate  the  prospects  of  monolithieally  implementing  a  nondispersive 
delay  line  on  a  given  YKS-film  strip  by  changing  the  magnitude  and  direction  of  the 
internal  bias  field  within  the  delay  line  between  the  two  microstrip  transducers  so  that  a 
portion  of  the  delay  line  would  support  MSFVW  and  the  remainder  of  the  delay  line 
would  support  MSBVW.  This  configuration  constitutes  a  monolithic  implementation  of 
the  Sethares  et  al  geometry  [l]  wherein  a  discrete  forward  wave-type  delay  line  was  cas¬ 
caded  to  a  discrete  backward  wave-type  delay  line,  with  the  requisite  dispersion  self¬ 
compensation  achieved  by  matching  the  delay  characteristics  of  the  two  delay  lines. 

A  detailed  study  of  MSFVW-to-MSBVW  mode  conversion  through  a  region  of 
bias-field  discontinuity  is  described  in  Chapter  1.  The  principal  objective  of  this  study 
was  to  establish  if  efficient  mode  conversion  is  possible  over  a  wide  bandwidth.  The 
problem  was  attacked  by  performing  theoretical  and  experimental  studies  simultane¬ 
ously.  The  experimental  study  indicated  the  presence  of  mode  conversion  but  this  con¬ 
clusion  was  somewhat  ambiguous  because  of  the  rather  unsophisticated  nature  of  the 
biasing  magnet  circuit  employed  and  the  consequent  lack  of  definition  of  the  profile  of 
the  internal  bias  field  along  the  length  of  the  delay  line.  Initial  theoretical  results  for 
mode  conversion  across  a  single  boundary  separating  two  regions  with  fixed  bias-field 
angles  were  encouraging  in  that  they  showed  that,  for  a  small  change  of  tilt  angle  across 
the  boundary,  efficient  mode  conversion  is  possible  over  a  wide  band. 

The  theory  of  mode  conversion  was  successfully  extended  to  a  practical  geometry 
with  a  gradual  bias-field  discontinuity  region.  This  geometry  was  modeled  as  a 
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segmented  or  layered  geometry  along  the  length  of  the  discontinuity  region,  i.e.,  it  was 
approximated  by  a  number  of  tilted  bias-field  layers  or  regions,  with  the  bias-field  orien¬ 
tation  or  tilt  angle  being  constant  within  each  region  and  varying  slowly  from  one  region 
to  the  next.  While  this  work  was  being  performed,  a  sophisticated  magnet  system  pro¬ 
viding  precision  bias-field  change  over  the  length  of  a  YIG  film  sample  was  being 
designed.  At  this  time,  the  theoretical  computations  on  mode  conversion  through  the 
layered  structure  ran  into  a  snag  which  was  ultimately  traced  to  the  existence  of  the 
unique  MSW  cut-off  condition.  A  careful  look  at  the  implications  of  this  cut-off  condition 
led  to  the  catastrophic  conclusion  that  mode  conversion  from  MSFVW  to  MSBVW  in 
going  through  a  bias-field  discontinuity  region  is  theoretically  not  possible. 

The  foregoing  turn  of  events  led  to  a  termination  of  the  effort  on  MSFVW  to 
MSBVW  mode  conversion  in  going  through  a  region  of  bias-field  discontinuity  region. 
The  effort  during  the  remainder  of  the  contract  period  was  dedicated  to  a  characteriza¬ 
tion  of  the  dispersion  of  an  MSFVW  pulse  suffered  during  propagation.  This  study  was 
successfully  performed  using  Fast  Fourier  Transform  techniques  and  is  described  in 
Chapter  2. 
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CHAPTER  1 

MSFVW-to-MSBVW  MODE  CONVERSION  AT  A  REGION  OF  BIAS  FIELD 

DISCONTINUITY 

1.1.  INTRODUCTION 

The  present  chapter  describes  theoretical  and  experimental  work  performed  on  the 
subject  of  MSFVW  to  MSBVW  mode  conversion  that  one  expects  to  occur  when  an 
MSFVW  is  incident  at  a  region  of  bias-field  discontinuity  separating  a  region  with  bias 
field  normal  to  the  YIG  film  (the  MSFVW  region)  from  a  region  with  bias  field  parallel 
to  the  YIG  film  and  to  the  direction  of  MSW  propagation  (the  MSBVW  region).  This 
study  was  undertaken  to  establish  that  the  expected  efficient  mode  conversion  is  indeed 
possible  over  a  wide  bandwidth  with  such  a  scheme.  The  aim  of  the  contract  was  to 
proceed,  after  establishing  the  existence  of  efficient  MSFVW-to-MSBVW  mode  conver¬ 
sion  through  a  suitably  profiled  bias-field  discontinuity  region,  to  a  monolithic  implemen¬ 
tation  of  a  nondispersive  MSW  delay  line  wherein,  within  a  single  delay  line,  part  of  the 
delay  line  supports  MSFVW  propagation  and  the  remaining  delay  line  supports  MSBVW 
propagation,  thereby  permitting  a  mechanism  for  achieving  dispersion  self-compensation. 
The  latter  geometry  constitutes  a  monolithic  implementation  of  the  Sethares  et  al  [l] 
nonditpenive  MSW  delay  line  geometry  which  consisted  of  two  cascaded  discrete  delay 
lines,  one  being  a  forward  wave-type  (MSSW  or  MSFVW)  delay  line  and  the  other  a 
backward  wave-type  (MSBVW)  delay  line. 

After  establishing  the  existence  of  efficient  mode  conversion  over  a  wide  bandwidth, 
the  question  of  achieving  dispersion  self-compensation  may  then  be  handled  in  the 
manner  it  was  handled  in  the  Sethares  et  al  cascaded  geometry  of  discrete  delay  lines, 
i.e.,  through  the  use  of  appropriately  spaced  ground  planes  (and  the  possible  use  of 
multi-layered  YIG  films  (2|).  Some  theoretical  work  performed  under  partial  support 


from  the  present  contract  (2]  (see  Attachment  A)  has  shown  that  superior  nondispersive 
MSFVW  characteristics  may  be  achieved  by  using  a  layered-YIG-film  in  conjunction 
with  a  spaced  ground  plane.  The  dispersion  self-compensation  is  not  treated  in  the 
present  report. 

The  theoretical  work  evolved  over  three  phases.  In  the  first  phase,  a  theory  for 
MSFVW-to-MSBVW  mode  conversion  was  worked  out  for  the  case  of  an  abrupt  bias- 
field  discontinuity.  The  choice  of  abrupt  bias-field  discontinuity  was  based  on  the  fact 
that  this  represents  the  simplest  geometry  and  the  consequent  expectation  that  this 
choice  might  yield  the  simplest  approach  at  handling  the  complex  problem  of  mode 
conversion.  While,  in  this  theory,  an  incident  fundamental-mode  MSFVW  was  con¬ 
sidered  to  scatter  into  all  MSFVW  reflected  modes  and  all  MSBVW  transmitted  modes, 
it  was  shown  that  the  resulting  infinite  set  of  equations  in  the  unknown  amplitude 
coefficients  could  be  truncated  depending  on  the  level  of  accuracy  desired  in  evaluating 
the  coupling  coefficients.  It  was  found,  in  agreement  with  what  one  might  expect,  that 
the  modes  that  need  to  be  included  in  the  truncation  are  the  ones  that  are  phase- 
matched  to  the  chosen  incident  mode  along  the  direction  normal  to  the  YIG  film  at  the 
boundary  separating  the  MSFVW  and  MSBVW  regions.  It  was  recognized  that,  in  addi¬ 
tion  to  the  multi-moded  nature  of  mode  conversion,  most  of  the  incident  power  is 
reflected  back,  a  feature  that  reflects  the  abrupt  nature  of  impedance  discontinuity  at 
the  junction  plane. 

Based  on  the  work  performed  during  the  first  theoretical  phase,  it  became  evident 
that  efficient  mode  conversion  into  a  single  mode  required  that  the  bias-field  discon¬ 
tinuity  be  sufficiently  gradual  so  that  the  transverse  phase-matching  to  the  incident 
mode  is  possible  with  essentially  a  single  converted  mode. 


An  extended  gradual  bias-field  discontinuity  region  may  be  approximated  by  a  seg¬ 
mented  or  layered  geometry  comprised  of  a  number  of  tilted  bias-field  regions,  with  the 
bias-field  orientation  or  tilt  angle  being  constant  within  each  region  and  varying  slowly 
from  one  region  to  the  next.  In  a  systematic  approach  at  solving  the  problem  of  mode 
conversion  through  such  a  “layered"  transition  region,  one  would  first  solve  the  problem 
of  mode  conversion  across  a  single  boundary  separating  two  regions  with  fixed  bias-field 
tilt  angles,  with  the  change  of  tilt  angle  across  the  boundary  being  small. 

The  problem  of  mode  conversion  across  a  single  boundary  was  successfully  solved 
for  the  case  of  an  incident  magnetostatic  volume  wave  (MSVW),  i.e.,  for  the  case  of  bias 
field  oriented  in  the  sagittal  plane  [3]  (see  Attachment  B).  It  was  found  that,  the  more 
gradual  the  discontinuity,  the  better  is  the  accuracy  in  using  a  three-mode  coupling 
model  which  assumes  the  coupling  of  three  modes,  viz.,  the  incident  wave,  one  reflected 
wave  and  one  transmitted  wave.  The  computations  have  shown  that  efficient  mode 
conversion  obtains  over  a  wide  band  for  modes  that  are  phase-matched  to  the  incident 
wave  in  the  direction  normal  to  the  YIG  film. 

In  the  third  theoretical  phase,  the  three-mode  coupling  theory  was  successfully 
extended  to  the  case  of  a  “layered”  transition  region  sandwiched  between  two  pure-mode 
regions  using  a  scattering-parameter  formalism.  In  the  course  of  numerical  computations, 
however,  irreconcilable  difficulties  were  encountered  which  were  ultimately  traced  to  the 
unique  cut-off  phenomenon  present  in  MSVV  propagation.  A  careful  evaluation  of  the 
implications  of  this  cut-off  condition  led  to  the  catastrophic  conclusion  that  mode 
conversion  from  MSFVW  to  MSBVW  in  going  through  an  extended  transition  region  is 
theoretically  not  possible. 

The  experimental  portion  of  the  work  was  focused  on  the  observation  of  mode 
conversion  in  a  single  straight  YIG  strip.  The  experiment  employed  a  3-transducer 
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geometry,  with  the  middle  transducer  placed  half-way  between  the  other  two  transduc¬ 
ers  being  a  diagnostic  transducer.  First  the  device  was  set  up  as  an  MSSW  delay  line 
within  an  electromagnet  and  the  transmission  was  measured  between  the  two  outer 
transducers.  The  choice  of  MSSW  for  the  forward  wave-type  portion  of  the  cascaded 
structure  was  based  on  the  fact  that  a  smaller  bias  field  is  required  compared  to 
MSFVW,  thereby  permitting  the  use  of  small  ceramic  magnets  in  the  experiment  and 
also  reducing  the  interaction  between  the  two  magnet  systems.  The  device  was  then  phy¬ 
sically  withdrawn  approximately  half-way  out  of  the  bias  field  so  that  a  transmitted  sig¬ 
nal  was  observed  at  the  center  transducer  but  no  signal  was  observed  at  the  outer 
receiving  transducer.  Next  a  horse-shoe  magnet  was  placed  over  the  withdrawn  region  so 
that  the  bias  field  orientation  for  this  region  was  for  MSBVW.  A  signal  was  now 
observed  at  the  outer  receiving  transducer  assumed  to  represent  the  converted  signal  at 
the  bias-field  discontinuity  region. 

The  main  problem  in  the  experiment  is  the  strong  interaction  between  the  two 
magnet  systems  because  of  their  proximity,  i.e.,  the  introduction  of  the  second  magnet 
produces  a  strong  perturbation  of  the  field  due  to  the  first  magnet.  A  precision  magnet 
design  providing  desired  bias  fields,  in  magnitude  and  direction,  in  the  two  pure-mode 
regions  and  separated  by  a  gradual  discontinuity  region  was  in  progress  at  the  time  of 
the  catastrophic  theoretical  discovery  that  MSFVW  to  MSBVW  mode  conversion 
through  a  region  of  bias- field  discontinuity  is  not  possible. 

While  the  design  of  a  precision  magnet  system  was  awaited,  a  special  magnet  cir¬ 
cuit  was  empirically  designed  and  machined  as  an  attachment  to  the  electromagnet  pole 
pieces.  The  objective  of  this  magnet  system  was  to  remove  the  problem  associated  with 
the  rather  large  dimensions  of  the  horse-shoe  pole  pieces  relative  to  the  length  of  the 
YIG  film. 
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A  jig  was  designed  and  machined  which  permitted  the  position  of  the  center  diag¬ 
nostic  transducer  to  be  variable  over  most  of  the  length  of  the  YIG  film.  An  important 
planned  diagnostic  study  utilizing  this  diagnostic  transducer  was  the  measurement  of  the 
physical  lengths  of  the  forward  and  backward  wave  regions  as  well  as  the  length  of  the 
transition  region. 


1.2.  THEORY 

The  abrupt  bias-field  geometry  (Fig.  1)  is  a  special  case  of  the  more  realistic  gra¬ 
dual  bias-field  geometry  (Fig.  2).  As  a  first  step  in  solving  the  complicated  problem  of 
mode  conversion  in  the  geometry  of  Fig.  2,  a  much  more  basic  problem  that  would  allow 
a  building-block  approach  was  solved.  The  two-region  geometry  shown  in  Fig.  3  is 
treated  first  wherein  the  internal  bias  fields  in  the  two  regions  are  oriented  in  the  sagittal 
plane.  The  building  block  approach  allows  the  results  for  the  geometry  of  Fig.  3  to  yield 
the  mode-conversion  properties  of  a  gradual  bias-field  discontinuity  region  since  the 
latter  may  be  approximated  by  a  number  of  tilted-bias- field  regions,  with  the  tilt  angle 
changing  slowly  from  one  region  to  the  next  as  shown  in  Fig.  2.  A  synthesis  procedure 
based  on  a  scattering-parameter  approach  has  been  employed  in  generating  the  results 
for  the  gradual  bias-field  discontinuity  geometry  from  the  single-interface  geometry. 

As  a  first  step  in  solving  the  single-interface  geometry  problem,  the  properties  and 
field  expressions  of  the  free  MSVW  modes  of  an  infinite  YIG  film  are  derived.  Then,  in 
solving  the  single- interface  geometry  problem  (Fig.  3),  which  assumes  a  single  MSWV 
mode  incident  at  the  boundary  separating  the  two  regions  with  different  bias  fields  in 
the  sagittal  plane,  one  postulates  the  reflected  and  transmitted  signals  to  be  given  by  a 
superposition  of  an  infinite  number  of  eigenmodes  of  the  two  respective  regions.  The 
unknown  amplitude  coefficients  of  these  eigenmodes  are  determined  by  the  mode- 
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Fig.  1.  Abrupt  bias-field  discontinuity  geometry 
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matching  technique,  i.e.,  by  applying  boundary  conditions  at  the  interface  plane  and 
invoking  the  hi  -orthogonality  [4,5]  properties  of  the  eigenmodes. 


1.2.1.  FREE  MSVW  MODES  OF  A  Y1G  FILM  WITH  INTERNAL  BIAS 
FIELD  LYING  IN  AN  ARBITRARY  DIRECTION  IN  THE  SAGITTAL 
PLANE 

The  YIG-film  geometry  with  a  tilted  internal  bias  field  lying  in  ihe  sagittal  plane  is 
shown  in  Fig.  4.  The  field  expressions  for  the  magnetic  potential  functions  representing 
the  free  MSW  modes  of  this  infinite  film  geometry  are  given  by  the  solutions  of  the  mag¬ 
netostatic  wave  equations 


which  satisfy  the  boundary  conditions,  i.e.,  the  continuity  of  tangential  magnetic  field 
and  normal  magnetic  induction  across  the  interfaces  :»0  and  The  potential  func¬ 

tions  *ajr,+yig  and  *ccc  representing  the  air,  YIG  and  GGG  regions,  respectively, 
yield  the  wave  rf  magnetic  fields  in  the  respective  regions  through  the  relation 
h=*-;ra<f  i>.  The  small-signal  permeability  tensor  JJ  for  the  rase  of  a  YIG  medium  biased 
beyond  saturation  by  an  internal  field  H,  lying  at  an  angle  8  to  the  r-axis  may  be 
readily  shown  to  be  given  by 
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Fig.  4.  Infinite  YIG-fibn  geometry  with  arbitrarily 
oriented  bins  field  ia  the  sagittal  plane 
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The  quantities  /  .  /,— 7<i, /f,  and  fy=iv,  M,  are,  respectively,  the  signal,  gyro,  and 
magnetization  frequencies. 

For  straight-crested  waves  propagating  along  s  direction,  i.e.,  d/dy=0,  the  poten¬ 
tial  functions  in  the  three  regions  have  the  expressions 
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and  the  propagation  wave  number  is  the  solution  of  the  dispersion  relation 


The  quantity  t  has  the  value  »— 1  for  *> 0  and  #  »  1  for  k  <0.  The  harmonic  time 
dependence  e  ' is  assumed  throughout  for  all  ac  quantities. 

The  transcendental  nature  of  the  dispersion  relation  in  Eq.  (6)  yields  an  infinite 
number  of  modal  solutions 
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where  m  is  a  modal  index  assuming  any  integer  value,  i.e.,  m  =0,1, 2.3,....  The  dispersion 
relation  (6)  and  the  field  expressions  (4)  reduce  correctly  to  the  more  commonly  used 
cases  corresponding  to  0-*  0'  and  00* . 

The  qualitative  features  of  the  MSVW  dispersion  diagram  corresponding  to  an  arbi¬ 
trary  angle  8  and  a  given  mode  are  shown  in  Fig.  5.  The  overall  passband  is  the  fre¬ 
quency  range  /,  </  <!  i  Within  this  passband,  the  wave  is  of  the  forward-wave  type 
for  frequencies  above  / and  of  the  backward  wave-type  for  frequencies  below  / The 
effect  of  changing  the  angle  8  on  the  dispersion  curve  for  a  given  mode  is  shown  in  Fig. 
6.  The  dispersion  curve  correctly  reduces  to  that  for  the  pure  MSFVW  obtaining  for 
J-0'  and  to  that  for  the  pure  MSBYW  obtaining  for  >■(»' . 

Notice  that  the  negative-*  branch  is  chosen,  as  it  must,  to  represent  the 
backward-wave  solution  rather  than  the  usually  shown  positive-*  branch.  The 
negative-*  branch  corresponds  to  a  positive  group  velocity  and  negative  phase  velocity 
and  is  the  appropriate  one  to  choose  for  the  geometry  of  Fig.  3  where  the  converted 
MSBYW  must  carry  energy  away  in  the  positive-z  direction.  (While  the  positive-* 
branch  for  MSBYW  is  a  valid  solution  of  the  dispersion  solution,  it  represents  energy 
propagation  in  the  negative-z  direction  and  thus  does  not  apply  in  the  present  problem.) 


The  total  time-average  power  carried  by  the  mth  mode  in  the  z  direction  is  given 


Fig.  S.  Qualitative  features  of  the  MSVW  dispersion  diagram 
for  aa  arbitrary  bias>leld  orieatatioa  aagie  aad  aay  one  mode 
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Fig.  0.  Elect  of  changing  «  on  the  MSVW  d»pe»ion  diagram  fot  a  given  mode. 


in  watts  per  unit  length  along  y  by  the  expression 
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where 

PB=iRe(e1,xh«*)  (9) 

is  the  Poynting  vector  for  the  mth  mode.  The  electric  field  e*  appearing  in  Eq.  (9)  is 
obtained  from  the  Maxwell  curl*  equation,  i.e., 

curt  *—-jun,n  h  (10) 

The  orthogonality  properties  of  the  guided  MSVW  modes  may  be  derived  based  on 

a  general  biorthogonality  relation  for  a  guiding  structure  containing  anisotropic  materials 
[4,5].  For  the  present  MSVW  modes,  this  biorthogonality  relation  translates  into  the  fol¬ 
lowing  four  relations: 
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The  superscripts  +  and  •  represent  modal  propagation  in  the  positive  and  negative  z 


directions,  respectively.  The  notation  used  to  represent  subscripts  is  that  a  quantity 
such  as  e„  represents  the  z  component  of  the  vector  ea ,  i.e.,  the  electric  field  associated 
with  the  mth  mode. 


1.2.2.  MSVW  MODE  CONVERSION  IN  A  SINGLE-INTERFACE 
GEOMETRY 

1. 2.2.1.  PROBLEM  FORMULATION 

The  geometry  of  the  problem  treated  here  is  that  shown  in  Fig.  3  where  x=0 
represents  the  interface  plane  and  the  incident  MSVW  mode  propagates  from  left  to 
right  along  the  positive  x  -axis.  In  the  mode-matching  technique,  the  scattered  fields  on 
the  two  sides  of  the  interface  plane  are  constructed  by  superimposing  all  the  free  modes 
of  the  YIG  film  situated  to  the  respective  side  of  the  interface  plane.  Thus,  the  total 
transverse  fields  in  the  region  x  <0  have  the  expressions 
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and  in  the  region  z  >0  the  expressions 
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where  the  superscripts  i.  r  and  t  represent  the  incident,  reflected  and  transmitted 
waves,  respectively,  and  the  primed  and  unprimed  quantities  represent  the  modal  fields 
in  the  x  >0  and  z  <0  regions,  respectively.  The  incident  wave  is  assumed  to  be  the  pth 
MSVW  mode  of  the  guiding  structure  in  the  z  <0  region.  The  quantities  R%t  and  T,f 
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are  the  unknown  reflection  and  transmission  amplitude  coefficients  corresponding  to  the 
nth  mode  constituent  of  the  reflected  and  transmitted  waves,  respectively.  The  boun¬ 
dary  conditions,  i.e.,  continuity  of  tangential  electric  and  magnetic  fields  at  the  interface 
plane  x  =0,  yield  the  infinite  set  of  coupled  equations 

+  E  R*f  »  E  T„  <v7  ;  *  =0  (20) 

*i+E^*»=E^V;  ; *=o  (21) 

where  the  unknown  coefficients  are  found  by  employing  the  biorthogonality  conditions 
(11)  to  (14). 

By  multiplying  Eq*.  (20)  and  (21)  by  (&»)'  and  (e'yr«)*>  respectively,  integrating 
the  resulting  equations  with  respect  to  z  over  the  entire  z  axis  and  then  adding  these 
equations  one  obtains 
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The  quantity  (C|r'm  n  may  be  considered  as  the  coupling  coefficient  between  the  nth- 
mode  component  of  the  reflected  wave  and  the  mth-mode  component  of  the  transmitted 
wave.  Similarly  the  quantity  (C) r ^  represents  the  coupling  coefficient  of  the  incident 
pth  mode  to  the  transmitted  mth  mode.  Since  the  index  m  runs  over  an  infinite  number 
of  integer  values,  i.e.,  m>0,l,2,.,  one  has  an  infinite  number  of  simultaneous  equations 


20 


in  an  infinite  number  of  unknowns  Rmp  (n  *=0,1,2,...).  The  solution  to  these  equations 
may  be  obtained  by  truncating  the  set  of  equations  at  a  finite  number  N  so  that  N 
amplitude  reflection  coefficients  R„p ,  n  =0,1,2,..., N ,  exhibit  a  coupling  between  one 
another.  The  corresponding  transmission  coefficients  Tnp,  n  =0,1,2, ...,N ,  may  then  be 
obtained  from  known  R'np s  by  using  the  relation 


where 
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The  power  conservation  condition  is  expressed  as 
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where  R£p—  — and  —7-  are  the  power  reflection  and  transmission 

Pip  Psp 

coefficients,  respectively,  corresponding  to  a  pth  incident  mode  and  mth  free-mode  com¬ 
ponents  in  the  reflected  and  transmitted  waves.  The  power  conservation  condition  is 
used  to  select  a  proper  truncation  numbei  N  as  well  as  a  pointer  to  indicate  the  accu¬ 
racy  of  the  solution. 

1. 2.2.2.  COMPUTED  MODE  CONVERSION  CHARACTERISTICS 

While  in  a  rigorous  theory  one  must  consider  the  reflected  and  transmitted  waves 


to  be  comprised  of  an  infinite  number  of  appropriate  free  plate  modes,  in  practice  one  is, 
of  course,  interested  in  mode  conversion  to  a  single  mode  and  not  to  multiple  modes. 
Computations  show,  in  agreement  with  expectation,  that  mode  conversion  to  a  single 
mode  improves  with  improvement  in  phase  matching  in  the  direction  normal  to  the  YIG 
film  between  the  incident  MSVW  and  the  desired  transmitted  MSVVV  mode.  This 
phase-matching  condition  implies  that,  the  more  gradual  the  discontinuity,  the  better  is 
the  accuracy  in  using  a  three-mode-coupling  model  corresponding  to  .V  =0  which 
assumes  a  coupling  of  three  modes,  viz.,  the  incident  wave,  one  reflected  wave  and  one 
transmitted  wave.  The  three- mode-coupling  computations  presented  here  are  for  an 
incident  fundamental-mode  wave  (p  =0,  x<0),  a  reflected  fundamental-mode  wave 
(m  =0,  x  <0)  and  a  transmitted  fundamental-mode  wave  (g  =0.j  >0). 

Three  frequency  bands  may  be  qualitatively  identified  in  interpreting  the  frequency 
variation  of  the  mode-conversion  process.  In  Fig.  7,  MSVW  dispersion  diagrams  are 
shown  for  two  neighbouring  value  of  9.  One  immediate  observation  from  this  figure  is 
that,  while  in  frequency  bands  A  and  C  mode  conversion  will  occur  between  modes  of 
the  same  type  (backward  wave-type  in  band  A  and  forward  wave-type  in  band  C),  in 
band  B  mode  conversion  will  occur  between  an  MSFVW  and  an  MSBVW.  Phase¬ 
matching  is  consequently  expected  to  be  better  in  bands  A  and  C  than  in  band  B.  Thus, 
one  may  expect  efficient  single-mode  conversion  over  a  wide  band  for  values  of  8 
approaching  0’  or  90".  On  the  other  end,  one  might  expect  mode  matching  problems 
for  frequencies  lying  within  band  B  or  for  frequencies  lying  within  bands  A  and  C  which 
are  close  to  the  resonance  frequencies  where  the  divergence  in  propagation  wave  number 
k  becomes  large.  These  predictions  are  basically  borne  out  by  the  computations. 

Computed  frequency  variation  of  the  MSVW  amplitude  reflection  R  and  transmis¬ 
sion  T  coefficients  is  presented  in  Figs.  8  to  12  for  different  values  of  0,  and  &2.  These 
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and  10*  of  6, 


mm 


computations  arc  for  a  YiG  film  of  thickness  20(im  and  an  internal  bias  field  magnitude 
of  1  kC  .  In  all  of  these  figures,  the  R  and  T  curves  are  shown  by  broken  and  dotted 
lines,  respectively. 

In  Fig.  8.  the  incident  MSWV  is  a  pure  MSFVW.  i  e  ,  0,—t)'  and  two  values  of 
are  taken,  viz.,  ^=5',  10'.  This  figure  shows,  in  agreement  with  expectation,  that  R 
goes  up  as  the  tilt-angle  discontinuity  is  increased  In  Figs  9  to  12.  the  value  of  0,  is 
increased  from  one  figure  to  the  next  so  that  the  entire  span  of  values  is  examined.  In 
these  figures,  the  R  characteristics  exhibit  a  characteristic  dip  at  a  frequency  of  about 
3.8  GHz.  It  is  apparent  that  phase-matching  must  be  good  at  this  frequency  Although 
it  is  not  obvious  why  this  should  happen  at  about  me  same  frequency  for  a  given  magni¬ 
tude  of  the  internal  bias  field  for  all  values  of  0,  and  02 .  computations  show  that  chang¬ 
ing  the  magnitude  of  the  internal  bias  field  in  either  region  may  relocate  the  frequency  of 
this  dip. 

In  the  figures  where,  within  the  B  band,  the  value  of  T  exceeds  unity,  an  attempt 
was  made  3t  improving  upon  the  three-mode-coupling  model  by  taking  into  the  coupling 
formalism  additional  modes.  The  computations  have  not  yielded  uniform  results  and  it 
would  generally  appear  that  the  series  expansion  technique,  i.e.,  the  mode-matching  tech¬ 
nique,  becomes  invalid  at  these  frequencies. 

The  computations  of  Figs.  8  to  12  are  devoid  of  irregularities  present  in  the 
corresponding  computations  of  Ref.  3.  The  latter  irregularities  manifest  themselves  m 
the  form  of  unacceptably  large  values  of  R  and  T  obtaining  at  the  high  and  low  fre¬ 
quency  bounds  of  the  MSVW  spectrum  and  possibly  at  other  frequencies  aho  It  turns 
out  that  the  difficulties  in  Ref  3  arose  from  the  use  of  orthogonality  conditions  that  are 
not  as  accurate  as  the  ones  used  in  this  report,  i.e  these  difficulties  are  not  attributed  to 
a  failure  of  the  three- mode-coupling  model. 


The  convergence  property  of  the  solution*  relative  to  the  number  of  modes  taken 
into  account  in  the  coupling  formalism  has  been  explored.  Sample  results  are  shown  in 
Tables  I  to  3  where  M  represents  the  number  of  modes  taken  into  account  in  the  i  >0 
region,  and  only  the  reflection  and  transmission  coefficients  corresponding  to  conversion 
into  fundamental  modes  are  given.  The  computations  for  each  table  are  for  fixed  values 
of  9X  and  #2  and  three  values  of  /lying  within  bands  A,  B  and  C.  The  rate  of  conver¬ 
gence.  in  general  agreement  with  expectation,  is  rapid  for  medium  values  of  the  propaga¬ 
tion  wave  number  k  and  mode  conversion  between  modes  of  the  same  wave-type. 

The  improvement  in  mode-matching  with  increase  in  the  number  of  modes  taken 
into  account  in  the  coupling  formalism  is  also  evident  in  the  computations  presented  in 
Fig.  13  where  modal  amplitudes  are  plotted  as  a  function  of  mode  number  for  a  201- 
mode-coupling  model.  The  modal  amplitudes,  in  agreement  with  expectation,  tend 
toward  smaller  values  as  the  modal  number  is  increased. 

1.2.3.  MSVW  SCATTERING  AT  A  REGION  OF  GRADUAL  BIAS-FIELD 
DISCONTINUITY 

The  gradual-bias-field-discontinuily  transition  region  is  modeled  as  a  cascade  of  S 
uniform  subregions  with  a  total  of  A'  +  l  interfaces  as  shown  in  Fig.  14.  In  the  ith  subre¬ 
gion,  which  extends  over  a  length  l, ,  the  bias-field  tilt  angle  has  a  uniform  value  of  9,  at 
all  points.  The  total  length  of  the  transition  region  is 

N 

i  «4 

1. 2.3.1.  SCATTERING  MATRIX  FOR  A  SINGLE  SUB-REGION  BASED 
ON  THREE-MODE-COUPLING  FORMALISM 

The  scattering  characterisation  of  the  itk  sub-region  (see  Fig.  15)  may  be  cast  in 
terms  of  three  scattering  matrices,  two  of  which  represent  the  scattering  properties  of 
the  boundary  planes  p  (i  )  and  ?(*  +  !)  of  the  sub-region  and  the  third  represents  the 
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Table  1  Illustrating  the  convergence  property  of  scattering  solution  with  the  number  of  modes 
taken  into  the  coupling  formalism  M  stands  for  the  number  of  reflected  modes  taken  into  account 
which  is  the  same  as  the  number  of  transmitted  modes  taken  into  account  The  computations  are 
for  t#, —  30f  and  $2=40' ,  i  e  ,  band  B  lies  between  3  36 GHz  and  3  67 GHz  . 


f 


3  3  GHz 
MSBVVV 

MSBVW 


(  Pul.  +  P,....)/  P," 


3  6  GHz 
MSF\  W 

MSBVW 


|  3  9  GHz 

1 

01943 

99802 

99643 

'  MSFVW 

2 

02039 

99796 

99927 

{ 

V 

1 

3 

02097 

99797 

1  00063 

• 

j 

i  MSFVW 

| 

4 

02138 

99799 

1  00159 

i 

1 

5 

02170 

99801 

1  00233 

F 

/ 

* 

i 

6 

02197 

99803 

1  00295 

> 

i 

! 

1 

02219 

99806 

1  00349 

A 

1 

8 

02238 

99808 

1  00396 

Table  2  Illustrating  the  convergence  property  of  Mattering  solution  with  the  number  of  modes 
taken  into  the  coupling  formalism.  M  stands  for  the  number  of  reflected  modes  taken  into  account 
which  is  the  same  as  the  number  of  transmitted  modes  taken  into  account.  The  computations  are 
for  and  50* ,  i.e.,  band  B  lies  between  3.83(7//*  and  3.09(7//* . 
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Fig.  13.  Modal  amplitude  reflection  coefficients  Rs0 ,  N=1  to  100,  plotted  as  a  function  of  modal 
index  N 
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property  that  the  length  of  the  sub-region  may  be  treated  as  a  transmission  line.  The 
amplitudes  of  the  incident  'waves  at  the  sub-region  are  taken  to  be  a,  and  a,+,  and  the 
amplitudes  of  the  reflected  waves  to  be  4,  and  41+I,  with  the  latter  taking  into  account 
the  mode  conversion  that  occurs  between  the  planes  p  (» )  and  p  (i  +1). 


A.  SCATTERING  MATRIX  FOR  A  SINGLE  INTERFACE 

The  scattering  process  at  a  single  interface  p(n)  separating  two-infinite  regions  n 
and  m  is  shown  in  Pig.  16.  The  signal  a ,  incident  from  region  n  scatters  into  a  reflected 
wave  Rt  an  and  a  transmitted  wave  Tmn  a,.  Similarly  the  signal  am  incident  at  the  inter¬ 
face  from  region  m  scatters  into  a  reflected  wave  Rm  o„  and  a  transmitted  wave  Tnm  am  . 
Thus  a  total  of  four  waves  are  involved,  i.e.,  the  incident  waves  a„  and  a„  and  the 
scattered  waves  6„  and  bn .  These  four  waves  are  related  by  the  scattering  matrix 
expression 


where 


8(»  i — 


Su  5,2 


The  subscript  n  of  the  scattering  matrix  refers  to  the  boundary  plane  p(n).  If  am  =0, 
the  problem  is  reduced  to  the  case  of  three-mode  coupling.  The  reflection  and  transmis¬ 
sion  coefficients  R„ ,  Rm,  Tnm  and  Tmn  were  obtained  in  Section  1.2.2.  Thus,  S(„)  is  com¬ 
pletely  known  for  any  plane  p  (n  )  based  on  the  solution  of  the  three-mode-coupling  prob¬ 
lem.  The  general  case  which  also  includes  higher-order  modes  has  also  been  treated  (see 
Section  1. 2.3.3.). 

B.  SCATTERING  MATRIX  FOR  TRANSMISSION  LINE 

The  input  and  output  quantities  of  a  transmission  line  representing  the  length  ln  of 
sub-region  n  are  shown  in  Fig.  17.  From  elementary  transmission-line  theory  these  quan¬ 
tities  are  readily  seen  to  be  related  by  the  expression 


(32) 


where 


S((„)= 


(33) 


The  quantities  k  ±  are  the  wave  numbers  associated  with  the  waves  propagating  in  the 
±z  -directions. 


C.  SCATTERING  MATRIX  FOR  A  SINGLE  SUB-REGION 


The  scattering  characterization  of  a  single  sub-region  is  thus  given  by  the  cascade 


relationship 


(34) 

where  S(1+I),  S<t)  and  S,(()  are  the  scattering  matrices  representing  the  planes  p(t+l)  and 
p  (« )  and  the  guide  length  l,  of  the  sub-region,  respectively. 

1. 2.3.2.  SCATTERING-MATRIX  FORMALISM  FOR  THE  TRANSITION 
REGION  BASED  ON  THREE-MODE-COUPLING  MODEL 

The  recognition  that  the  overall  transition  region  represents  a  cascade  of  individual 
sub-regions  yields  the  overall  transition-region  scattering-matrix  relationship 


«.+i 


S(.+i)SMl)S(l) 


aw+> 

[bN+l 


=  W 


(35) 


where 
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In  Eq.  (35),  ax  is  the  incident  MSFVW  amplitude,  b ,  is  the  reflected  MSFVW 
amplitude  and  6w+l  is  the  transmitted  MSBVW  amplitude.  If  it  is  assumed  that  the 
MSBVW  region  is  a  semi-infinite  one  or  matched  at  the  end,  then  o/v+l=0.  The  overall 
transition-region  amplitude  reflection  and  transmission  coefficients  may  be  related  to  the 
components  of  W ,  i.e., 
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1.2. 3. 3.  GENERALIZED  MULTI-MODE  COUPLING  FORMALISM  FOR 
THE  TRANSITION  REGION 

Consider  a  single-interface  plane  p(n)  which  separates  two  sub-regions  n  and  m 
where  m=n+l  (see  Fig.  18).  The  incoming  waves  at  the  interface  have  amplitudes  a, 
and  am  and  the  outgoing  waves  have  the  amplitudes  o,  and  bm  .  These  amplitude  quan¬ 
tities  are  now  defined  by  the  column  matrices 


o"=  (a„  0  .  a„ ,  , 
om=  («™o  ■  i  , 
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•  .  «.N) 
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where  the  superscript  T  represents  the  operation  of  transpose.  Each  element  in  the 
column  matrices  has  two  subscripts,  the  first  one  of  which  denotes  the  region  and  the 
second  one  the  mode  number.  Thus,  am  is  the  ith  incident  mode  in  region  n  and  bmj  is 
the  jth  outgoing  mode  in  region  m.  The  a 's  and  f  s  are  related  by  the  reflection  and 
transmission  matrices  given  by 
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The  elements  of  the  matrices  in  Eq.  (41)  are  interpreted  as  follows.  /?,"  represents 
the  reflection  coefficient  corresponding  to  the  plane  p  (n )  which  is  the  ratio  of  the  itk 
outgoing  mode  amplitude  to  the  jtk  incident  mode,  both  modes  being  in  the  nth  region. 
Similarly,  T,'m  is  the  transmission  ratio  of  the  itk  outgoing  mode  on  the  n  side  to  the  jtk 
incident  mode  on  the  m  side.  All  these  elements  are  determined  in  a  similar  analysis  to 
that  in  Section  1.2.2.  The  scattering  properties  of  the  interface  plane  p(n)  are  given  by 
the  expression 
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is  a  2(,V  +  l)x2(N  +  1)  matrix.  On  each  side  of  an  interface  plane,  a  total  of  N+l  incident 
and  N+l  reflected  modes  are  assumed. 

The  transmission  property  of  a  section  of  the  guide  of  length  /,  (see  Fig.  19)  may 
be  represented  as 
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where 
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where  t,+  and  k.  are  the  wave  numbers  for  the  tth  positively  propagating  and  the  )th 
negatively  propagating  modes  respectively,  in  the  /,  guide  section. 

The  overall  mode  conversion  property  of  the  entire  transition  region  may  be 
expressed  as 
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The  W  matrix,  which  is  a  2(.V+l)x2(A/  +  l)  matrix,  fully  determines  the  mode 
conversion  properties  of  the  transition  region,  i.e.,  the  overall  transition-region  reflection 
and  transmission  coefficients  can  be  computed  from  a  knowledge  of  the  incident  mode  or 
modes  in  the  pure  MSFVW  region  Thus,  if  it  is  assumed  that  the  incident  wave  is  a  sin¬ 
gle  mode  which  is  the  fundamental  MSFVW  mode  and  the  MSBVW  region  is  semi¬ 
infinite  (or  matched),  then 
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where  «l0  is  the  amplitude  of  the  fundamental  MSFVW  mode.  The  reflection  coefficient 
to  the  ilk  MSFVW  mode  due  to  the  entire  transition  line  is 
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The  transmission  coefficient  to  the  jth  MSBVW  mode  is 
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The  expressions  for  R,0  and  T,0  can  be  readily  put  in  terms  of  the  elements  of  the  W 
matrix.  Numerical  computations  of  mode  conversion  through  the  transition  region  may 
then  be  performed. 


1.3.  EXPERIMENTAL  WORK 

The  experimental  work  evolved  over  three  phases.  In  the  first  phase.  MSW  delay¬ 
line  fabrication  techniques  were  perfected  resulting  in  delay  lines  with  a  strong  suppres¬ 
sion  of  direct  pick-up  (typically  in  excess  of  50  dB  over  most  of  the  2  GHi  to  12  GHi 
band)  and  a  low  insertion  loss  (typically,  without  tuning,  in  the  order  of  5  dB  in  the 
MSSW  configuration  and  10-15  dB  in  the  MSFVW  or  MSBVW  configuration). 

The  delay  lines  employed  microstrip  transducers  defined  on  Cu-clad  Duroid  soft 
substrates  which  were  found  to  work  as  well  as  those  on  Au-clad  alumina  hard  sub¬ 
strates  and  are  yet  much  cheaper.  YlG-film  strips  with  ends  beveled  with  diamond  pol¬ 
ishing  were  used  and  a  high  level  of  ripple  suppression  in  the  pass  band  was  achieved, 
particularly  for  the  case  of  MSSW’s. 
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Id  the  secoDd  experimental  phase,  a  basic  experiment  on  mode  conversion  was  per¬ 
formed.  A  3-transducer  geometry  was  employed,  with  the  middle  transducer  placed 
half-way  between  the  other  two  transducers  being  a  diagnostic  transducer.  The  device 
was  first  set  up  as  an  MSSW  delay  line  and  the  transmission  was  measured  between  the 
two  outer  transducers.  The  device  was  then  physically  withdrawn  approximately  half¬ 
way  out  of  the  bias  field  so  that  a  transmitted  signal  was  observed  at  the  center  trans¬ 
ducer  but  no  signal  was  observed  at  the  outer  receiving  transducer.  The  observed  ampli¬ 
tude  and  phase  characteristics  of  the  signal  picked  up  at  the  center  transducer  are  shown 
in  Fig.  20.  Next  a  horse-shoe  magnet  was  placed  over  the  withdrawn  region  so  that  the 
bias-field  orientation  in  this  region  was  for  MSBVW.  A  signal  was  now  observed  at  the 
outer  transducer  (see  Fig.  21).  The  delay  characteristic  seen  in  Fig.  21  corresponds  to  a 
BYW  and  is  in  qualitative  agreement  with  the  overall  delay  characteristic  of  a  cascade  of 
an  MSSW  and  MSBVW  delay  lines  at  the  low-frequency  end  of  the  overlap  frequency 
band 

f  or  the  cascaded  geometry  comprised  of  an  MSSW  delay  line  followed  by  a 
MSBVW  delay  line,  one  actually  expects,  for  the  case  of  overlapping  passbands,  the 
delay  characteristic  to  initially  show  a  decrease  in  delay  time  with  increasing  frequency 
and  then  an  increasing  delay  time  with  further  increase  in  frequency.  The  observed 
delay  characteristic  in  the  monolithic  structure  exhibits  a  basically  monotonic  decrease  in 
delay  time  with  increasing  frequency,  a  feature  that  is  at  variance  with  expectation  at 
the  high  frequency  end  of  the  passband  spectrum.  It  was  concluded  that,  while  the 
observed  signal  in  Fig  21  might  represent  mode  conversion,  this  was  not  necessarily  fully 
proved. 

The  difficulties  in  the  experiment  lie  mainly  in  the  implementation  of  a  controlled 
slowly  varying  bias  field  in  the  transition  region  as  well  as  the  inevitably  strong  interac- 


Fig  20.  Amplitude  and  phase  characteristics  observed  at  the  center  transducer  when  only  the 
MSSW  bias  field  is  present  and  the  delay  line  ts  withdrawn  half-way  out  of  the  magnet  poles. 


tion  between  the  magnet  systems  because  of  their  close  proximity.  Furthermore,  the  size 
of  the  poles  of  the  horse-shoe  magnet  is  quite  significant  compared  to  the  dimensions  of 
the  YIG  film.  Thus,  while  it  seemed  like  a  good  first  try  to  withdraw  the  YIG  film  half¬ 
way  out  of  an  MSSW  bias  field  and  then  impose  the  field  of  a  horse-shoe  magnet,  the 
need  for  a  careful  design  of  a  precision  single  magnetic  circuit  was  fully  recognized. 

In  order  to  perform  a  more  careful  study  of  the  mode-conversion  phenomenon  in 
the  above  experiment,  a  three-transducer  delay-line  device  was  constructed  which  per¬ 
mitted  the  position  of  the  center  diagnostic  transducer  to  be  variable  over  most  of  the 
length  of  the  YIG  film.  A  photograph  of  this  device  is  shown  in  Fig.  22.  An  important 
diagnostic  study  that  was  attempted  albeit  without  success  is  that  of  measuring  the  phy¬ 
sical  lengths  of  the  forward  and  backward-wave  regions  as  well  as  the  length  of  the  tran¬ 
sition  region.  The  difficulty  in  performing  this  study  arose  from  the  sensitivity  of  the 
two  interacting  magnet  systems  to  any  slight  movement  of  either  magnet  system.  One 
would  draw  the  conclusion,  as  a  result  of  the  second  phase  of  the  experimental  study, 
that,  while  MSW  transmission  in  a  delay  line  subjected  to  two  biasing  magnetic  systems 
(one  intended  to  produce  an  MSSW  bias-field  orientation  in  the  first  half  of  the  delay 
line  and  the  other  intended  to  produce  an  MSBVW  bias-field  orientation  in  the  second 
half  of  the  delay  line)  was  observed,  it  was  not  conclusive  that  unambiguous  mode 
conversion  had  indeed  been  observed. 

The  third  phase  of  the  experimental  study  was  concerned  with  the  design  and 
implementation  of  a  precision  magnetic  circuit  which  was  machined  out  of  steel  and 
screwed  on  to  the  electromagnet  poles.  While  such  a  precision  magnet  design  was 
awaited,  an  empirical  design  was  machined  and  attached  to  the  poles  of  the  electromag¬ 
net.  A  photograph  of  this  empirical  design  is  shown  in  Fig.  23.  A  photograph  of  the 
new  magnet  poles  with  the  three-transducer  delay  device  placed  within  them  is  shown  in 


Fig.  24.  The  bias-field  orientation  profile  was  tested  using  iron  filings  and  was  qualita¬ 
tively  of  the  type  desired.  A  tuning  of  the  orientation  or  magnitude  of  the  bias  field  was 
achieved  with  the  use  of  coils  placed  around  the  poles.  Experimental  observations  were 
made  which  were  inconclusive.  However,  before  further  work  could  be  carried  out,  new 
theoretical  understanding  into  the  mode-conversion  process  was  attained  which  indicated 
a  serious  problem  in  the  achievement  of  mode  conversion  (see  section  1.4)  and  led  to  a 
termination  of  the  general  effort. 

An  alternative  approach  to  the  problem  of  observing  mode  conversion  had  also 
been  initiated  which  involved  the  use  of  a  curved  YlG-film  strip,  e.g.,  a  right-angle  bend. 
The  main  attraction  of  using  a  curved  YIG-film  strip  is  that  it  should  afford  a  much 
better  control  over  the  bias  fields  separating  the  two  wave-type  regions.  A  central  ques¬ 
tion  attaching  to  the  practicality  of  using  a  YIG-film  bend  whether  the  propagation  loss 
in  going  around  the  bend  is  tolerable.  The  phosphoric-acid  process  for  etching  out  a 
curved  YIG-film  strip  has  been  well  established  in  our  laboratory  but  the  theoretical 
difficulties  explained  in  Section  1.4  also  forced  a  termination  of  this  effort. 

1.4.  NEW  THEORETICAL  DEVELOPMENTS 

Computations  of  mode  conversion  in  a  tingle-interface  geometry  with  slightly 
differing  bias  fields  on  either  side  of  the  interface  have  shown  efficient  mode  conversion 
over  a  wide  band  (see  Section  1.2.2.).  However,  attempts  at  calculating  the  mode- 
conversion  characteristics  for  the  case  of  a  gradual  bias- field- discontinuity  transition 
region  ran  into  snags  which  were  tracked  down  to  the  k=0  point  on  the  MSVW  d.sper- 
sion  diagram  (see  Fig.  6)  occuring  at  a  unique  value  of  6  for  a  fixed  frequency.  The  k=0 
point  represents  a  cut-off  condition:  A  lower  cutoff  frequency  exists  even  when  the  elec¬ 
tromagnetic  retardation  effect  is  taken  into  account  (see  Attachment  C).  At  the 


Fig.  24.  Photograph  of  added  magnet  poles  with  the  3-transdncer  delay-line  device  placed  within 
them. 


point  along  the  guide  where  cut-off  obtains  propagation  along  the  guide  simply  ceases  to 
exist.  In  other  words,  the  sub-region  of  the  transition  region  where  the  bias-field  tilt 
angle  0  for  the  ehoten  frequency  yields  the  MSW  wave  number  1=0  simply  prevents  the 
signal  from  going  through  it.  It  follows  that  mode  conversion  from  MSFVW  to  MSBVW 
in  going  through  a  gradual  bias-field-discontinuity  transition  region  is  theoretically  not 
possible.  This  catastrophic  condition  has  been  investigated  in  considerable  detail  in  an 
attempt  at  finding  possible  means  of  overcoming  the  cut-off  condition  but  this  effort  has 
been  unsuccessful. 

1.5.  CONCLUSIONS 

A  theoretical  formalism  for  handling  the  problem  of  MSVW  mode  conversion  across 
a  region  of  bias-field  discontinuity  was  successfully  completed  employing  the  mode¬ 
matching  technique.  Computations  of  mode  conversion  across  a  plane  of  small  discon¬ 
tinuous  change  in  bias-field  orientation  indicate  efficient  mode  conversion  over  a  wide 
bandwidth.  Attempts  at  computing  the  mode-conversion  properties  of  a  gradual-bias- 
field-discontinuity  region  ran  into  a  snag  which  was  tracked  down  to  the  existence  of  the 
1=0  point  on  the  MSVW  dispersion  diagram.  For  any  value  of  the  signal  frequency 
within  the  MSV'W  passband,  there  exists  a  unique  angle  lying  between  0*  and  90*  which 
yields  1=0  for  the  MSVW  wave  number.  The  1=0  point  represents  a  catastrophic  condi¬ 
tion  since  it  implies  that  MSVW  propagation  and,  therefore,  conversion  would  not  occur 
across  the  sub-region  of  the  transition  region  where  1=0.  This  catastrophic  discovery 
rendered  further  work  on  MSVW  mode  conversion  across  a  region  of  bias-field  discon¬ 
tinuity  without  merit  and  thus  forced  a  termination  of  the  effort. 


CHAPTER  2 


DISPERSION  OF  A  TIME-LIMITED  MSFVW  CW  PULSE 
2.1.  INTRODUCTION 

MSW’s  are  pronouncedly  dispersive  in  geometries  devoid  of  dispersion-control 
features.  The  present  chapter  reports  the  main  results  of  a  detailed  theoretical  and 
experimental  study  of  the  dispersion  of  a  time-limited  MSFVW  CW  pulse.  The  theory 
has  incorporated  the  frequency  response  of  the  input  and  output  microstrip  transducers 
in  3  delay-line  geometry  and  has  yielded  results  that  are  in  excellent  agreement  with  the 
experimentally  observed  oscilloscope  pictures  of  delayed  pulses  picked  up  by  an  output 
microstrip  transducer  whose  spacing  from  the  input  transducer  is  made  variable  through 
the  placement  of  the  output  transducer  on  a  platform  movable  with  a  screw  arrange¬ 
ment  (see  Fig.  1.22). 

The  theoretical  computations,  if  performed  by  a  brute-force  method,  would  be 
based  on  the  standard  Fourier  transform  technique  Thus,  by  representing  a  time- 
limited  CW  pulse  in  terms  of  its  Fourier  spectrum  one  simply  finds  the  change  in  the 
amplitude  (if  losses  are  included)  and  the  phase  of  every  single  frequency  component  of 
the  signal,  and  then  at  the  end  of  the  propagation  the  individual  frequency  components 
are  combined  to  synthesize  the  delayed  pulse  signal 

Difficulties  in  the  inverse  transform  process  arise  whenever  the  integrand  of  the 
transformation  is  an  ill-behaved  function,  e  g.,  a  rapidly  oscillating  function  which  is 
very  much  the  case  with  MSW  pulses  An  elegant  approach  to  overcoming  this  problem 
is  provided  by  the  fast  Fourier  transform  (FFT)  technique  which  assumes  that,  instead 
of  a  solitary  or  single  pulse,  there  exists  a  train  of  pulses  that  are  periodic  in  time.  The 
existence  of  a  train  of  pulses  sufficiently  separated  in  time  is  incidentally  always  the  case 
in  experiments, 
i 
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The  FFT  technique,  with  its  attendant  process  of  summing  over  weighted  discrete 
Fourier  components,  makes  the  computational  process  extremely  efficient,  i.e.,  the  com* 
putational  time  is  significantly  shorter  because  of  fewer  computational  points  required 
than  in  the  brute-force  or  direct  method  of  performing  the  integration  represented  by 
the  inverse  Fourier  transform.  The  degradation  of  the  accuracy  in  the  computations 
goes  up  as  the  number  of  computational  points  per  temporal  cycle  goes  down.  This 
degradation  of  envelope  distortion  is  not  large  with  the  FFT  technique  even  when  only 
one  computational  point  per  carrier  cycle  is  employed.  The  FFT  technique  also  allows 
propagation  loss  to  be  included  in  the  calculation  since  complex  quantities  are  an  intrin¬ 
sic  part  of  the  FFT  approach.  A  detailed  theory  of  the  FFT  technique  as  well  as  its 
implementation,  described  in  a  1985  doctoral  thesis  at  Stony  Brook  [6],  is  not  presented 
here. 

2.2.  BASIC  THEORY  OF  PULSE  DISPERSION 


Let  the  CW  pulse  be  defined  by  the  function 

=0,/ )  =  .4  n(-^)  cos  2*/c  t  (1) 

at  its  starting  point  x=0  where  A  is  the  pulse  amplitude,  / c  is  the  carrier  frequency  and 

n(y)  is  the  window  function 


niit-l1'  1  <TI  2 

r  |0,  t  >r/2 


(Strictly,  a  precisely  rectangular  MSFVW  input  pulse  is  not  possible  because  of  the  finite 
bandwidth  over  which  MSFVW’s  exist). 


Since  the  function  $z=0,«)  is  symmetric  in  time  and  extends  only  over  the  time 


period  |  t  |  <r/2,  the  Fourier  transform 


of  c>(j  =0,f )  may  be  written  as 


*(*  =0 ,/  )  =  (/  -ft )r+«ne  (/  +/,  )rj  (4) 

Since  NlSFVWs  exist  only  in  the  band  f,<f  </  s,  only  positive  frequencies  are  allowed 
in  Eq.  (4),  i.e., 

*(*=0./  )  =  —sine  (f  -fc  )t  (5) 

The  foregoing  signal  <j>(z=0,t)  was  at  the  point  x=0.  As  time  increases,  the  pulse  will 
move  so  that  at  any  point  x>0, 


$(*■/  )=  *(z=0J  )exp  (  jkz)  (6) 

The  inverse  Fourier  transform  of  Eq.  (6)  reconstitutes  the  pulse  in  space  and  time,  yield¬ 
ing 


<t<(z  ,t  )  =  A  T 


' r  sin  *(/  -fc)r 

fK  (f-fc)r 
1  » 


cos  (2jt  / 1  -kz)  df 


(7) 


2.3.  NUMERICAL  RESULTS 

The  computational  efficiency  of  the  FFT  technique  vis-a-vis  the  direct  or  brute- 
force  integration  technique  in  the  computation  of  pulse  distortion  is  demonstrated  in  the 
computations  of  Figs.  1  and  2.  The  FFT  computations  are  given  in  Fig.  1  while  the 
direct  integration  results  are  given  in  Fig.  2.  These  figures  represent  the  shape  of  a  20 
ns  pulse  of  rf  frequency  / c  =4  GHz  after  propagation  through  a  distance  of  1  cm.  Three 
sets  of  computations  are  presented  in  each  of  these  figures  corresponding  to  2,4  and  8 
computed  points  per  carrier  cycle,  respectively.  The  shape  of  the  pulse  envelope  in  Fig. 
1  remains  essentially  unchanged  as  the  number  of  computed  points  per  cycle  is  reduced 
from  8  to  2.  On  the  other  hand,  the  computations  employing  the  direct  integration 
approach  (which  computes  the  carrier  signal  rather  than  the  envelope)  show  a  significant 
degradation  of  the  envelope  as  the  number  of  computed  points  per  rf  cycle  is  reduced. 


5' 


Fig.  1.  FFT-computed  envelop*  of  a  20  ns  MSFVW  pulse  of  rf  frequency  f  c  -=4  GHz  after  a  pro¬ 
pagation  pathlength  of  1cm  for  three  different  values  of  computed  points  Nc  per  rf  cycle:  (a) 

Nt  -2,  (b)  Ne  -4,  (c)  Nt  -8. 


so 


in  Fig.  3,  the  spreading  of  the  envelope  of  a  pulse  of  width  20  ns  at  x=0  is  shown 
for  different  positions  x  of  the  observation  point,  i.e.,  x=0.5  cm,  0.75  cm,  1.00  cm  and 
1.25  cm.  As  expected,  the  distortion  occurs  primarily  in  the  form  of  stretching  of  the 
pulse  width  and  the  attendant  “peaking"  of  the  pulse. 

A  significantly  smaller  level  of  envelope  distortion  achievable  in  a  nondispcrtivc 
geometry  is  demonstrated  in  the  computations  of  Fig.  4  which  are  for  a  nondispersive 
double-YlG-film-layered  structure.  In  the  computations,  the  upper  YIG  film  has  a  thick¬ 
ness  of  65  pm  and  saturation  magnetization  of  2kG,  the  lower  film  has  a  thickness  of  10 
pm  and  saturation  magnetization  of  1.75  kG,  and  a  ground  plane  is  placed  at  a  distance 
of  100pm  above  the  top  YIG  film.  The  overshoot  at  the  skirts  of  the  pulse  envelope  is 
the  unavoidable  result  of  the  finite  bandwidth  over  which  MSFVWs  exist. 

2.4.  EXPERIMENTAL  RESULTS 

The  experimental  work  on  the  dispersion  of  an  MSFVW  pulse  was  carried  out  using 
a  41pm  thick  YIG  film  with  a  saturation  magnetization  of  about  1.5kG.  The  delay  line 
measurements  made  at  rf  frequencies  /«  *=2.06  GHZ  and  /c=2.16  GHz  are  shown  in 
Figs.  5  and  6,  respectively.  These  measurements  were  made  directly  on  the  Tektronix 
7104  oscilloscope  without  the  use  of  a  crystal  detector.  It  is  evident  that,  except  for 
spurious  signals  following  the  main  delayed  pulse,  the  agreement  between  experiment 
and  theory  is  generally  excellent. 

2.5.  CONCLUSION 

The  results  of  a  theoretical  and  experimental  study  of  MSFVW  pulse  distortion 
with  propagation  are  reported.  It  is  demonstrated  that  the  FFT  technique  of  computa¬ 
tion  of  the  distortion  in  pulse  envelope  with  propagation  is  significantly  more  efficient 


1.25cm 


J. 


i 


0.75cm 


J 


i.5cm 


than  the  direct  or  brute-force  integration  technique.  The  agreement  between  experiment 
and  theory  on  the  shape  of  the  delayed  pulse  is  generally  excellent. 
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FI*.  1  ■  Geometry  of  the  problem 


Nondispenive  Magnet  ostatic-Forward-Voiuine-Wave 
Propagation  in  a  Double-Layered  YIG-Fibn  Structure 

).  P.  PAREKH  and  K.  W.  CHANG 

A  detract- Computed  group -delay  chancteriftics  an  preaeated  lO ul¬ 
tra  tin*  the  significant  improvement  in  bandwidth  for  aondMpnnive 
magne lot tabc-forward-volume -wave  propaption  that  may  be  achieved 
utihring  a  douMe-YlG-film-layercd  structure  over  a  singie-YIG-film 
geometry. 

A  theoretical  prediction,  reported  briefly  elsewhere  |l|,  that  a 
significant  improvement  in  bandwidth  for  nondispershe  magnetostatic- 
forward-volume-wave  (MSFVW)  propagation  over  the  single-YlG-fQm/ 
dielectric/ground-plane  configuration  may  be  achieved  using  a  double- 
YIG-film -layered  structure  in  conjunction  with  a  separated  ground  plane 
is  described  here  tor  saturation  magnetization  parameters  that  are  readily 
realizable  in  practice.  The  computed  delay  versus  frequency  character¬ 
istics  presented  here  furthermore  show  that  (he  requirement  for  the 
achievement  of  bandwidth  improvement  stated  in  (1|,  Le..  the  YIG 
film  closer  to  the  ground  plane  possess  the  higher  saturation  magnetiza¬ 
tion  of  the  two  films,  may  in  fact  be  relaxed  so  as  to  permit  either 
film  to  possess  the  higher  saturation  magnetization. 

The  geometry  of  the  problem  under  consideration  is  shown  in  Fig.  1 . 
The  YIG  film  ot  thickness  dy  grown  on  the  GGG  substrate  (Y1G-1  film) 
is  overlaid  with  a  different  YIG  film  of  thickness  <fj  (Y1G-2  film).  A 
peitcct  electrical  ground  plane  is  placed  at  a  height  h  above  the  top  film 
surface.  (Although  a  practical  ground  plane  would  give  rise  to  conduc¬ 
tive  loss,  the  computations  show  that  optimum  MSFVW  nondispersive 
characteristics  are  obtained  if  the  ground  plane  is  sufficiently  separated 
from  the  top  -YIG-film  surface  (typically  by  100  um  or  more)  so  that 
the  conductive  loss  is  rendered  insignificant  (in  the  order  uf  1  dB/cm  or 
less)).  An  external  saturating  bias  field  Ho  *  #o*o  applied  normal  to 
the  films  produces  internal  fields  Hu  =  H0  -  Mox  iQ  and  H^  mB0- 
•**03*o  ,n  the  two  films  where  Moi  and  Mo j  are  the  values  of  the  satu¬ 
ration  magnetization  for  the  respective  films.  For  straight-crested 
MSFVW  propagation  along  x,  with  x  and  r  variation  exp/(2*/r  -  kx), 
the  dispersion  relation  is  ( 1 1 

fur  ♦  mi  tanh  kh )  tan  0|*d|  tan  0j *dj  ♦  0((m2  ♦  tanh  kh)  tan  0* kdj 

*  0z(mi  ♦  tanh  kk)  tan  ffykdy  *  0,0a(l  ♦  tanh  kh)  ■  0  (I) 

where  u*  *  -0*.  with  k  *  1,  2.  is  the  xx  or  yy  component  of  the  rela¬ 
tive  permeability  tensor  characterizing  the  *th  YIG  film.  The  expression 
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the  U.S.  Air  Force  (Rom*  Air  Development  Center)  under  Contract 
FI9628-82-K-0009 

The  authors  are  with  the  Department  of  Electrical  Engineering,  State 
University  of  New  York  at  Stony  Brook,  Stony  Brook,  NY  11794. 
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Fig.  2.  Computed  variation  of  group  delay  with  frequency.  The  solid 
curves  represent  the  case  where  the  lower  saturation-magnetization 
film  (uqAfox,  3  1 300  G)  is  closer  to  the  ground  plane  and  the  dashed 
curves  represent  the  case  where  the  higher  saturation-magnetization 
film  ( Mo  Mqh  *  1750  G)  Is  closer  to  the  ground  plane  The  fixed  pa¬ 
rameters  in  the  computations  are  h  «  127  Mm.  (thickness  of  lower 
saturation-magnetization  film)  «  10  pm.  and  external  bias  field  mo  Ho  ' 
2500  G.  The  thickness  dy/  of  the  higher  saturation-magnetization  film 
is  used  a*  a  variable  parameter. 


for  mx  is 

M  k-(/J  -/!*)/(/’-/'?*)  <2> 

where  fok  *  7M0Hft  and  /j*  *  l/o*</o*  */*/*>! 1/2  *rt,  respectively  . 
the  lower  and  upper  frequency  bounds  of  ihe  MSFVW  spectrum  for  the 
individual  k th  film.  The  frequencies  fox  znd  f\rk  •  7Mo-**ok  ere.  re¬ 
spectively.  the  gyro-  and  magnetization  frequencies  for  the  *th  film 

The  dispersion  relation  (I)  correctly  reduces  to  that  for  MSFVW ’s  in 
the  single-YlG-film/dielectnc/ground-plane  structure  in  the  limit  of  zero 
thickness  of  either  film.  The  resulting  two  single-YlG-film  structures 
have  dispersion  curves  that  are  qualitatively  similar  but  shifted  in  fre¬ 
quency  By  applying  the  usual  mode  coupling  criteria  at  the  crossover 
points  of  these  two  sets  of  MSFVW  single-film  dispersion  curves  one 
readily  deduces  the  qualitative  features  of  the  MSFVW  dispersion  char¬ 
acteristics  for  the  double- YIG-film -layered  structure  which  are  in  agree¬ 
ment  with  computed  characteristics.  The  present  letter  highlights  only 
those  aspects  of  the  dispersion  characteristics  of  MSFVW's  in  a  double- 
layered  structure  which  show  that  a  significant  improvement  in  nondis- 
perstve  MSFVW  propagation  is  possible  utilizing  such  a  structure  rather 
than  a  single-YlG-fUm  geometry. 

The  computations  show  that  the  MSFVW  propagation  is  nondispersive 
over  a  frequency  band  located  at  the  lower  end  of  the  spectrum  over 
which  the  variation  of  MSFVW  fields  across  both  films  is  sinusoidal  in 
nature,  this  spectrum  being  defined  by  the  joint  condition  mi  <  0. 
mj  <  0.  This  implies  that  a  nondispersive  MSFVW  solution  k  to  disper¬ 
sion  relation  (I)  exists  over  a  band  of  frequencies  lying  immediately 
above  the  larger  of  the  two  gyro-frequencies /0(  and  foi  .  The  computa¬ 
tions  show  that  a  bandwidth  of  more  than  1 0  percent  with  a  delay  ripple 
of  less  than  tOJ  ns  may  be  readily  achieved  by  an  appropriate  choice 
of  parameters  d ( ,  dj ,  h.  Mot  •  Mot  •  Ho •  All  the  computations  pre- 
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FI*.  3.  Computed  group-delay  duractarltUc*  for  fixed  parameters 
b  =  92  tun.  dg,  *  10  «un,  HqHq  =  2500  G,  M gifgg,  m  1 500  G,  and 
U0iV0/f  *  1750  G,  with  dff  used  a*  a  variable  parameter.  Tha  dia- 
tlnction  betwaan  tha  aoiid  and  dashed  curraa  ia  tha  aama  aa  in  FI*.  2. 


tented  here  are  for  the  lowest  or  fundamental  MSFVW  mode  of  the 
structure.  The  values  of  the  parameters  used  in  the  computations  are 
indicated  in  the  captions  of  Figs.  2-4. 

The  improvement  in  MSFVW  nondispersive  propagation  that  results 
from  using  a  double-YlG-fUm-layeted  structure  is  apparent  in  Fig.  2 
where  the  frequency  variation  of  group  delay  is  presented,  with  the  thick¬ 
ness  of  the  higher  saturation-magnetization  film  taken  as  a  parameter. 
The  effect  on  the  group-delay  characteristics  of  the  saturation  mag¬ 
netization  of  the  YIG  film  closer  to  the  ground  plane  having  the  higher 
or  the  lower  value  of  the  two  films  is  also  shown  in  Fig.  2. 

The  difference  in  the  group-delay  characteristics  resulting  from 
whether  it  is  the  higher  or  the  lower  saturation-magnetization  film 
that  is  placed  closer  to  the  ground  plane  is  highlighted  on  an  expanded 
group  delay  scale  in  Figs.  3  and  4.  The  computations  of  Fig.  3  are 
quazi -optimized  to  yield  constant- delay  characteristics  for  the  case 
where  the  higher  saturation-magnetization  film  is  placed  closer  to  the 
ground  plane.  The  group-delay  characteristics  which  obtain  when,  for 
the  same  parameters,  the  locations  of  the  two  films  relative  to  the 
ground  plane  are  interchanged  are  also  shown  in  Fig.  3. 

A  choice  of  different  parameters  a  required  in  order  to  obtain  quazi- 
optimized  constant-delay  characteristics  with  the  lower  saturation- 
magnetization  film  placed  closet  to  the  ground  plane  (see  Fig.  4).  It 
is  evident  from  the  computations  of  Fgs.  3  and  4  that  the  optimum 
parameters  for  MSFVW  nondispersive  propagation  are  changed  if  the 
locations  of  the  two  films  relative  to  the  ground  plane  are  interchanged. 


2.  a  2.9  3.  a  3.  :  3.  2 


Fig.  4.  Computed  group-delay  chanctarietka  for  fixed  parameters 
h  *  127  urn,  dg,  *  10  jun,  pqHq  «  2500  G,  momOL  *  G>  *nd 
hqMoh  3  >750  G,  with  d/f  uaad  aa  a  variable  parameter.  The  distinc- 
tkw  between  the  solid  and  dashed  curves  is  the  same  aa  in  Fig.  2. 


In  conclusion,  the  present  computations  employing  readily  realizable 
values  of  saturation-magnetization  for  the  two  (Urns  have  shown  that 
MSFVW  nondispeishre  characteristics  similar  to  those  reported  in  ( 1  ] 
(where  p0M0  values  of  1750  and  2000  G  were  used)  may  be  achieved 
without  requiring  to  engage  in  the  difficulties  of  growing  a  low-loti 
epitaxial  YIG  film  with  hqM0  *  2000  G.  It  has  also  been  shown  that 
the  YIG  film  that  is  closer  to  the  ground  plane  may  have  either  the 
higher  or  the  lower  saturation  magnetization  of  the  two  films  but  opti¬ 
mum  parameters  (or  MSFVW  nondispersive  propagation  ate  not  main¬ 
tained  if  the  locations  of  the  two  films  relative  to  the  ground  plane  are 
interchanged. 
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ABSTKAl ! 

The  present  paper  reports  theoretical  computations 
of  node  conversion  which  arises  from  the  scatter¬ 
ing  of  a  magnetostatic  volume  wave  in  a  YIG  film, 
with  internal  bias  field  oriented  in  an  arbitrary 
direction  in  the  sagittal  plane,  at  a  boundary 
separating  two  regions  with  different  bias-field 
orientations,  a  three-mode-coupling  nodel  is  em- 
pleved  based  on  the  assumption  of  a  weak  bias-field 
discontinuity.  Efficient  mode  conversion  is  found 
to  obtain  for  modes  that  are  phase-matched  in  the 
direction  normal  to  the  YIG  film  to  the  incident 
wave.  The  bandwidth  for  mode  conversion  is  largest 
when  the  bias  field  is  either  normal  or  parallel  to 
the  .  IC  fi  In. 


L-.  INTRODUCTION 

t'nlike  SAW's,  magnetostatic  waves  (MSW's)  exhibit 
a  pronounced  dispersive  behaviour.  This  feature 
makes  MSW's  more  difficult  to  exploit  for  signal¬ 
processing  applications.  However,  MSW  dispersion 
control  techniques  do  exist  and  much  effort  has 
1  been  expended,  and  is  currently  directed,1 at 

|  the  realization  cf  viable  nondispersive  as  well  as 

,  linearly  dispersive  delay  lines. 

'  i  i.e  principal  MSW  dispersion  control  technique  con- 

i  sists  of  placing  a  ground  plane  in  proximity  to  the 

i  YIG  film,  thereby  changing  the  MSW  field  profile 

and  modifying  the  dispersion  characteristics.  An¬ 
other  promising  approach  at  dispersion  control  is 
the  use  oi  a  multiple- layered  YIG  film  rather  than 
\  t!ie  usual  single  YIG  film.  Theory  predicts  signi- 

,  f  leant  Improvement  in  dispersion  control  when  a 

|  lavered  YIG  film  is  used  in  conjunction  with  a 

|  suitablv  spaced  ground  plane. 

i 

1  Sethares  et  al*  have  demonstrated  that  a  nondis¬ 

persive  MSW  delav  line  may  be  realized  through 
dispersion  self-compensation  that  is  obtained  by 
cascading  two  linearly  dispersive  MSW  delay  lines, 
one  corresponding  to  a  forward  wave-type,  i.e., 
either  magnetostatic  surface  wave  (MSSW)  or  magne¬ 
tostatic  forward  volume  wave  (MSFVW) ,  and  the  other 
to  a  backward  wave-tvpe,  i.e.,  magnetostatic  back¬ 
ward  volume  wave  (MSBVW) .  The  realization  of  a 
nondispersive  delav  line  with  acceptable  pc-For»- 

*Work  supported  by  l.S.  Air  Force  (RADC)  under 
Contract  No.  K18628-82-0009 


n.y.  ir- 

ance  for  some  applications  has  recently  been  re¬ 
ported  which  employs  a  ground  plane  with  variable 
spacing  in  the  forward-wave-tvpe  (MSSW)  portion  of 
the  cascaded  structure  and  a  grouno  plane  with 
fixed  spacing  in  the  MSBVW  portion  of  the  struc¬ 
ture.  The  variable  spacing  of  the  ground  plane  ir. 
the  MSSW  delay  line  is  reported  to  provide  the 
needed  improvement  in  the  linearity  of  dispersion 
of  trie  delay  line. 

The  present  work  is  motivated  by  the  conjucture 
that  a  monolithic  implementation  of  the  Sethares 
et  ai  cascaded  geometry  might  be  possible  on  a 
single,  YIG  film.  The  monolithic  approach  requires 
mode  conversion  from  MSFVW  (say)  to  MSBVW  through 
the  use  of  different  bias-field  orientations  in  the 
two  VIG-film  regions  (see  Fig.  1).  The  problem  ad¬ 
dressed  in  the  present  paper  is  that  of  determining 
if  efficient  node  conversion  into  a  single  mode  is 
possible  across  a  bias-field  discontinuity.  The 
problem  of  establishing  the  level  of  dispersion 
self-compensation  is  the  same  as  that  applying  to 
the  Sethares  et  al  cascaded  geometry  and  is  not  ad¬ 
dressed  here. 

A  gradual  bias-field  discontinuity  is  more  realis¬ 
tic  in  practice  than  an  abrupt  bias-field  discon¬ 
tinuity  (see  Fig.  2).  It  turns  out  that  a  gradual 
bias-field  discontinuity  is  actually  necessary  for 
efficient  mode  conversion  into  a  single  mode.  A 
slow  change  in  bias-field  tilt-angle  and  magnitude 
ensures  that  an  incident  fundamental -mode  MSFVW  is 
phase -matched  essentially  to  a  single  converted 
mode,  i.e.,  the  fundamental-mode  MSBVW. 


YIG 

[  !>-* 
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MSf-vw  *)>•*  •  user*  ryio* 

Figure  1.  Abrupt  bias-field  discontinuity  geometry. 
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As  a  first  step  in  solving  the  complicated  problem 
of  mode  conversion  across  a  bias-field  discontin¬ 
uity  region,  a  much  more  basic  problem  which  uill 
allow  a  building-block  approach  has  been  solved. 

The  two-region  geometry  shown  in  Fig.  3  is  treated 
wherein  the  internal  bias  fields  in  the  two  regions 
are  oriented  in  the  sagittal  plane.  The  building 
block  approach  will  allow  the  present  results  to 
yield  the  mode-conversion  properties  of  a  gradual 
bias-field  discontinuity  region  since  the  lacter 
nay  be  approximated  by  a  number  of  tilted-bias- 
field  regions,  with  the  tilt  angle  changing  slowly 
from  one  region  to  the  next. 

In  a  rigorous  theory,  one  considers  an  incident 
magnet . static  volume  wave  (MSVW)  to  scatter  into 
al_l  reflected  MSVW’s  and  all  transmitted  MSYN's. 
this  means  that  there  will  he  an  infinite  set  of 
equations  in  the  unknown  amplitude  coefficients  of 
the  scattered  waves,  While  a  systematic  procedure 
for  the  evaluation  of  these  amplitude  coefficients 
ha^  been  established  which  is  based  on  an  ortho¬ 
gonality  condition  between  modes  of  a  given  wave- 
tvpe,  the  real  interest,  of  course,  is  in  mode 
conversion  to  a  single  mode  and  not  to  multiple 
modes . 


Computations  show,  in  agreement  with  expectation, 
that  node  conversion  to  a  single  mode  improves 
with  improvement  in  phase  matching  in  the  direction 
normal  to  the  YIG  film  between  the  incident  >!SW 
and  the  desired  transmitted  MSVU  mode.  This  phase- 
matching  condition  implies  that,  the  more  gradual 
the  discontinuity,  the  hotter  is  the  accuracy  in 
using  a  three-mode-coup  1 ing  model  which  assumes  a 
coupling  of  three  modes,  viz.,  the  incident  wave, 
cr.e  reflected  wave  and  one  transmitted  wave. 

I  he  present  work  is  confined  in  scope  to  the  three- 
rode-coupling  model.  A  weak  bias-field  discontin¬ 
uity  is  assumed,  with  the  bias  field  changing  by  a 
small  amount,  in  direction  only,  in  going  through 
the  discontinuity  plane.  Computed  amplitude  reflec¬ 
tion  and  transmission  charac teri st ics  are  presented 
which  demonstrate  the  general  validity  of  the  three- 
mode-coupling  model  provided  the  bias-field  discon- 
• i au i t v  i s  weak . 
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Figure  3.  The  node-conversion  geometrv  treated  in 
this  paper. 


II.  MSVW'S  IN  A  YIG  FILM  WITH  TILTED  INTERNAL 
BIAS  FIELD 


The  YIG-film  geometry  with  a  tilted  internal  bias 
field  lying  in  the  sagittal  plane  is  shown  in  Fig. 
4.  It  is  well-known  that  the  YIC  film  supports  a 
pure  MS F W  when  0=0°  and  a  pure  MSBVW  when  8»90’  . 

It  has  been  shown  by  Weinberg’  that,  for  an  inter¬ 
mediate  value  of  8,  the  MSVW  dispersion  diagram 
consists  of  two  branches,  one  corresponding  to  a 
forward  wave-type  and  the  other  to  a  backward  wave- 
type.  These  branches  meet  at  k=0  and  a  frequence 
determined  hv  the  value  of  0.  The  latter  frequency 
represents  a  transition  between  a  forward  wave-tvpe 
and  a  backward  wave-type,  with  the  upper  and  lower 
dispersion  branches  exhibiting  a  resonance  at  the 
upper  and  lower  bounds  of  the  MSVW  spectrum,  res¬ 
pectively. 

The  MSVW  dispersion  diagrams  for  two  neighbouring 
values  of  0  are  shown  in  Tig.  5.  It  is  seen  that, 
while  in  frequency  bands  A  and  C  mode  conversion 
will  occur  between  modes  of  the  same  tvpe  (forward 
wave-type  In  band  C  and  backward  wave-tvpe  in  band 
A),  in  band  B  mode  conversion  will  occur  between 
an  MSFVW  and  an  MSBVW.  Phase-matching  is  conse¬ 
quently  expected  to  be  good  in  bands  A  and  C  but 
poor  in  band  B.  Thus,  one  may  expect  efficient 
mode  conversion  over  a  wide  band  for  values  of  «1 
approaching  0°  or  90°  but  significant  difficulties 
in  the  validity  of  the  three-mode-coupling  model 
might  occur  for  frequencies  lying  within  band  B. 

The  latter  difficulties  do  not  matter  much  for 
values  of  0i  approaching  0°  or  90°.  However,  for 
values  of  approaching  45°,  the  frequency  band  B 
will  lie  somewhere  near  the  middle  of  the  MSVW 
spectrum.  One  might  also  expect  difficulties  in 
achieving  phase-matching  in  the  three-mode-coupling 
model  for  frequencies  approaching  the  resonance 
frequencies  where  the  divergence  in  propagation 
wave  number  k  becomes  large.  These  predictions  are 
basically  borne  out  by  the  computations. 


Figure  ?.  Gradual  bias-field  discontinuity  geometry. 
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Figure  5.  MSW«'  dispersion  curves  for  two  neigh¬ 
bouring  bias-field  angles. 


vs1",  dispersion  curves  for  an  arbitrar- 
1  lv  oriented  bias-field  in  the  sagittal 

p !  me . 


I  ,  ,'Wirf.f’  MODE  CONVERSION  CHARACTERISTICS 

.-nputed  frequency  variation  ot  the  MSVU  reflection 
g  and  transmission  T  coefficients  is  presented  in 
Figs.  6  to  !1  for  different  values  of  8^  and  *2- 
; hese  computations  are  for  a  V X 0  film  of  thickness 
3  )_m  and  an  internal  bias  field  magnitude  of  lkC. 

In  ilg.  h,  the  incident  MSV'W  is  a  pure  MSFVV,  i.e., 
and  8j  is  increased  from  2.5”  to  10°.  This 
figure  shows,  in  agreement  with  expectation,  that  R 
goes  up  as  the  angle  discontinuity  is  increased. 

The  difficulty  in  the  validity  of  the  three-mode 
coupling  model  at  the  high  and  low  frequency  bounds 
of  the  MSVK  spectrum  is  evident  in  the  computations 
where  the  R  value  becomes  unacceptably  large. 

In  Figs.  7  to  13,  the  value  of  is  increased  from 
one  figure  to  the  next  so  that  the  entire  span  of 
values  is  examined.  In  each  case,  the  value  of  02 
is  set  .lose  to  that  of  0j ,  with  the  differential 
Seo  rning  smaller  as  ©i  approaches  the  critical  angle 
ot  from  either  direction.  The  difficulty  asso- 

dated  with  hand  B  in  Fig.  5,  which  is  somewhat 
alleviated  bv  decreasing  the  differential  angle  as 
•*;  iprre.ic  hex  the  4  5*  angle,  is  reflected  In  Figs,  6 
t,.  Ll.  \  noticeable  effect  is  that  the  transmission 
bandwidth  decreases  as  01  is  increased  from  0*  tc 
-'i”  md  then  increases  as  01  is  further  increased 
*  r iS*  to  °nJ.  At  small  as  well  as  large  values 

f  •*[  ,  the  transmission  bandwidth  is  seen  to  be 
quite  large.  In  Figs.  6  to  13,  the  R  characterls- 
h, ,  exhibit  dips  at  frequencies  where  phase-match¬ 
ing  is  good. 
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Figure  6. 
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Frequency  variation  of  reflection  R  and 
transmission  T  coefficients  for  fixed 
0j“O"  and  variable  Si- 
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figure  “s.  Frequency  variation  >f  R  and  T  for  ftj -15® 
ami  ^2*^"* 


Frequence  variation  of  R  and  !  tor  ftj  =  30  ‘ 
and  ft 2 ~  30 •  3C  • 
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Figure  10.  Frequency  variation  of  R  and  T  for 
and  •  1  ° * 


1983  ULTRASONICS  SYMPOSIUM  -  241 


-requencv  variation  of  R  and  T  for  e.-bO' 
and  e-,-60.  3°  .  1 


f,  GH? 

Iigure  U.  frequency  variation  of  K  and  T  for  8j-75° 
and  **•>=  76°  . 
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! he  present  work  has  demonstrated  the  general  valid¬ 
ity  of  the  three-mode-coupling  model  provided  the 
bias-field  discontinuity  is  weak.  If  the  bias-field 
discontinuity  is  not  weak,  the  need  arises  to  take 
into  account  more  transmitted  and  reflected  modes. 
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Figure  13.  Frequence  variation  of  R  and  1  for 
and  02-90' . 

The  computations  have  shown  that  efficient  moot  con¬ 
version  obtains  for  modes  that  are  phase -ma t c bed  in 
the  direction  normal  to  the  YI<  film  to  the  incident 


The  bandwidth  for  mode  conversion  increases  si.-nfi- 
cantlv,  with  a  simultaneous  relaxation  in  the  need 
for  a  weak  bias-field  d  i  scont i nu  i  t y ,  as  the  an.  M 
api  roaches  0°  or  °0C  from  the  critical  ib°  value. 
Improved  bandwidth  for  values  of  f)j  approaching  .  i * 
requires  that  the  bias-field  discontinuity  is  ren¬ 
dered  sufficiently  weak.  The  frequency  band  h  ol 
Fig.  b  where  mode  mismatch  exists  and,  therefore, 
the  need  arises  to  take  into  account  more  modes  in 
the  oupling  formalism  than  dene  here  can  be  reduced 
in  bandwidth  through  the  choice  oi  a  suitablv  vu« 
bias-field  discontinuity  but  nit  full v  eliminated. 
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Dielectrically  induced  surface  waves  and  the  magnetodynamic 
modes  of  a  YIG  plate  7  ? 
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Hie  rni^.-ir i  \uri.wc  .ui.i  igiunic  ir<\!cs  o!  a  YKi  plate  metallized  at  one  surface  are  described  for 

the  case  of  propagation  at  right  -ngle*  to  a  ungentiaiiv  applied  bias  field  The  connection  of  the 
dtelecfncalli  .’rwjucexj  surface  »ave  on  a  Y  IG  halt-space  u>  the  magr.rtodv camic  -Plume  modes  of  a  YIG 
plate  is  uuluated 
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Parekh  and  Pan  inigi1  have  reported  the  propagation 
characteristics  of  a  raagr.ctodyn  tunic  surface  wave  which 
they  labeled  a  dielectrically  induced  surface  wave 
'DISW).  The  choice  >t  this  nomenclature  was  based  on 
the  fact  that  this  surface  wave  does  not  obtain  in  a 
magnetostatic  .walvs.s  awl  that  it  also  ceases  to  obtain 
it.  the  limit  -  i  where  is  the  relative  permittivity 
of  the  YIG  substrate.  In  this  communication,  we  wish 
to  pr;:it  out  tli.it  if.  instead  of  the  half-space  geometry 
used  in  Ref.  1  .  a  YIG  plate  of  finite  thickness  is  con¬ 
sidered,  the  dispersion  branch  for  the  DISW  becomes 
part  of  a  m<>di first  TF  dielectric  mode  dispersion  curve. 
Thus,  the  latter  nomenclature  may  actually  be  used  in 
place  of  the  nomenclature  DISW.  However,  since  the 
dielectric  modes  of  a  plate  disappear  in  the  limit  of  a 
half- -.pace,  the  appropriateness  of  the  nomenclature 
DISW  :s  not  invalidated.  This  distinction  car.  be  .•  nmifi- 
cant  even  though  in  practice  one  always  deals  with  a 
YIG  p'ate  rather  than  a  half-  -pace  if  a  resistive 
absori  ?r  is  placed  on  ..ne  surface  of  a  YTG  plate,  this 
c.  nf i.Cj ratioti  may  be  mor,  ar-mrc'eT.  modeled  1  >.  a 
YTG  half-space  thin  by  .1  YIG  plate. 

The  TE  magnctodynamic  nnides  of  a  tangc ntiallv  mag¬ 
netized  YIG  plate  have  recently  been  reported'-  *  for  the 
case  of  propagation  at  right  angles  to  the  litas  field. 
However,  the  study  of  '.jvr-hon  and  Nad  an1  (units  itself 
to  the  surface-wave  nodes,  and  thus  does  not  bring  out 
!ln  ltd  itu  tisfap  of  the  DISW  to  the  TF  magnetr-'.h  namic 
volume  mode,-,  of  a  YTG  piale  nor  to  the  dielectric 
modes  to  which  these  volume  modes  reduce  if  the  spins 
tn  YIG  are  considered  to  be  clamped.  Furthermore, 
the  ex.., -fence  of  a  bai  kward-  wave  characteristic  over 
par!  '  the  DISW  disper-inn  curve  pertiucr.t  to  the 
[>artfi!i\  e  ‘■tiUt.u’ri  YIG  plate  geometry,  whert-i  ’  one 
surface  of  the  YTG  plate  is  metallized  and  the  other 
adjoined  to  vacuum,  is  not  indicated  in  Ref.  3.  A  brief 
description  of  the  complete  spectrum  of  the  TE  mag¬ 
neto-1,  ..innr  modes  of  a  partial!-.-  metallized  YIG  plate 

is  if-.,  refore  Vi  e  h.-r  e 

f  ois1  i»r  ’he  Y  Hi  p.;ii  to  .crept,  the  space  C  i  i , 
with  in.  serf  ..re  i  i  oiai.-d  with  a  perfectly  c  vnt-  .ding 
metall.i  f>tn  ir-f  ’he  surface  t  it  adjoined  to  v  icu  im 
Assuming  tb.  I  i..s  field  f  .  be  oriented  along  the  .-  li_ 
rectum.  and  the  gi id«l  waves  to  propagate  along  v  with 
the  variation  exp>!2-r?  -  h  \ !,  the  dispersion  relation  for 
the  rmg  netrsfvnairr-  ’nodes  of  'he  YTG  plate  is  found  to 
be 

46.  Nu  M 


**„, I  ft-?  -  3  +  6*yl  tanlOtJ,  -  fc’)‘  '*  d |  s  (**,  -  =  0, 


(1) 


where  G  2m (ncfr)t/2f  is  the  wave  number  of  homo¬ 
geneous  plane  waves  in  vacuum,  and  k0/-  (fr/P,tf )‘ ,2k0 
is  the  wave  number  of  plane  waves  in  YIG  for  propaga¬ 
tion  at  right  angles  to  the  bias  field.  The  effective 
permeability  titll  ( P  -  f\ )  (  f2  -  t\)  is  the  reciprocal 
of  the  i»  or  iy  component  of  the  inverse  of  the  per¬ 
meability  tensor  characterizing  the  YIG  medium,  and 
jb  )/„  f  {f2  -  f\)  is  the  iv  (or  the  negative  of  the  v.i ) 
component  of  the  inverse  of  the  permeability  tensor. 
The  frequencies  f ,  and  /6  have  the  expressions  /, 

(/o(/o  *  'wM1  -  and  /,  f0  *  tu  .  with  fa  r_  •,  n9Ha  2r  being 
the  gyrof requency  and  -  )  p„.\f a  2ir  the  magnetization 
frequency 


The  solution  kv  of  Eq.  (1)  for  a  given  frequency  / 
corresponds  to  a  surface  wave  if  k2f  -  Jtj  <  0  and  to  a 
volume  wave  if  -  k*y  '  0.  Since  u,rr,  or  equivalently 
to,,  is  negative  within  the  frequency  band  /,</</,, 
only  surface  waves  can  obtain  within  this  band.  Outside 
of  this  band,  surface  waves  can  obtain  for  hf>k0f  and 
volume  waves  for  lt0  <  ky  <  .  Equation  (1)  reduces 

correctly  to  the  dispersion  relation  for  the  TE  dielectric 
modes  of  a  partially  metallized  dielectric  plate*  if  the 
spins  in  YTG  are  considered  to  be  clamped,  i.e.,  p#ff 
1  ario  b  0.  Also,  as  indicated  in  Ref.  3,  in  the  case 
of  surface  waves  Eq.  (1)  reduces  correctly  to  the  dis¬ 
persion  relation  for  magnetostatic  surface  waves  on  a 
partially  metallized  YTG  plate5  in  the  short-wavelength 
limit  t>v  k,:,  br,t.  and  to  the  dispersion  relation  for 
m.igtietodynamir  surface  waves  on  an  unrnplallized  YIG 
half  space1  3  in  the  limit  of  a  half-space,  i.e.,  d—  *• . 
Since  Eq  (1)  is  independent  of  whether  the  origin  of  the 
coordinate  system  is  on  the  metallized  or  the  unmetal- 
lized  surface  of  the  YTG  plate,  one  might  expect  that  in 
the  limit  of  a  half-space  the  surface-wave  case  of  Eq. 

<1!  should  decouple  into  two  dispersion  relations; 

:ir>ni“ly.  ihu  dispersion  relation  for  magnetodynamic 
surface  waves  on  an  unmetallized  half-space  x  •  0  and 
tin  dispersion  relation  for  magnetodvnamic  surface 
waves  on  a  metallized  half-space  r  *  0.  The  fact  that 
only  the  former  dispersion  relation  is  obtained  is  simply 
a  consequence  of  the  result  that  a  metallized  YIG  half¬ 
space  does  not  support  magnetodynamic  surface 
waves  * 


<  i  mputed  dispersion  diagrams  for  the  magnetody- 
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FIG.  1.  Computed  dispersion  diagram  (or  the  magnetodynamtc 
plate  modes  propagating  in  the  *y  direction  with  the  dc  bias 
field  In  the  ♦  *  direction.  The  solid  curves  represent  the  mag- 
netodynamlc  mode  solutions,  while  the  dashed  curves  corre¬ 
spond  to  the  TE  dielectric  mode  solutions  which  would  obtain  if 
the  spins  In  YIG  could  be  clamped.  The  dotted  curves  are  the 
dispersion  branches  for  plane  waves  in  YIG  progatlng  at  right 
angles  to  the  bias  field.  The  heavy  circles  Indicate  mode 
cutoff. 


namic  plate  modes  propagating  in  the  ±y  directions  are 
presented  in  Figs.  1  and  2,  respectively.  The  computa¬ 
tions  are  for  a  1 -cm -thick  YIG  plate  with  «rJ  =  16, 

-  1750  G  and  uo/f0  =  500  G.  For  the  purpose  of  compari¬ 
son,  the  TE  dielectric  mode  solutions4  for  the  same 
geometry  wherein  YIG  Is  assumed  to  be  magnetically 
clamped  are  also  shown  in  Fig.  1.  The  strong  effect  of 
the  spins  on  these  dielectric  modes  is  clearly  evident 

in  Figs.  1  and  2.  The  dispersion  curves  for  the  lowest 
three  magnetodynamic  volume  modes  corresponding  to 
the  values  of  the  mode  number4  m  =  0,  1 ,  and  2  are 
shown  in  Figs.  1  and  2.  For  each  value  of  m ,  where 
m  =  0,  1 ,  2,  . . . ,  two  dispersion  branches  for  the 
volume  modes  exist,  one  of  which  lies  above  /=/,  and 
the  other  below  /=/,.  The  frequency  /,  =  [f\ 

+  /„/«(!*  -  defines  the  crossing  point  of  the 

dispersion  curves  Af  =  A8  and  ky  -  A0/  The  volume 
modes  occurring  below  /=/,  exhibit  a  resonance  A, 

-  •  at  /=/,.  As  the  frequency  /  is  increased  to  /,  from 
below,  the  density  of  the  volume  modes  Increases 
mono  tonic  ally,  becoming  infinite  at  /=/,.  The  disper¬ 
sion  curves  of  all  the  volume  modes,  with  the  possible 
exception  of  the  dispersion  curve  of  the  upper  m  =  0 
mode  propagating  in  the  -  y  direction,  have  cutoff  fre¬ 
quencies  on  the  A?  =  A,  curve.  The  unique  features  of  the 
upper  m  =  0  volume  mode  propagating  in  the  -  y  direc¬ 


tion  will  become  apparent  following  a  description  of  the 
surface -wave  branches  in  Figs.  1  and  2. 

Only  one  surface-wave  branch  is  present  in  Fig.  1 
which  corresponds  to  the  modified  magnetostatic  sur¬ 
face  wave  “tied”  to  the  YIG-vacuum  interface.  The 
wave  number  ky  on  this  branch  increases  monotonically 
from  a  value  ky=  (/*//-),/*A#  at  the  lower  frequency 
bound  /,,  where  the  slope  df/dkt  (or  the  group  velocity) 
is  zero,  to  ky—  «  at  the  upper  frequency  bound  f4  =  /„ 

+  i/8.  In  contrast,  two  surface-wave  branches  are 
present  in  Fig.  2,  one  of  which  represents  the  modified 
magnetostatic  surface  wave  tied  to  the  metallized  YIG 
surface.  This  branch  is  tangent  to  the  dispersion  curve 
Af  =  A0  of  plane  waves  in  vacuum  at  the  lower  frequency 
bound  /,  (which  is  slightly  larger  than  /,  and  is  practical¬ 
ly  independent  of  the  plate  thickness),  and  exhibits  a 
resonance  ky~  «  ad  /=/,.  The  other  branch,  which 
represents  a  continuation  of  the  upper  m  =  0  magneto¬ 
dynamic  volume  mode  and  may  thus  be  identified  as  a 
modified  dielectric  mode,  lies  within  a  frequency  range 
ft  </</,,  where  ft  and  /,  are  both  functions  of  the  plate 
thickness.  In  the  limit  of  a  half-space,  ft  approaches 
/,  from  above,  and  ft  approaches  its  value  for  an  un¬ 
metallized  YIG  half-space1  from  below.  As  d  is  de¬ 
creased,  a  monotonic  increase  in  ft  and  a  monotonic 
decrease  in  /, take  place  until  for  a  critical  plate 
thickness  dcrtt  both  ft  and  fx  coalesce  at  a  frequency 


FIG.  2.  Computed  dispersion  diagram  far  the  mageetodynemlc 
plate  model  propagating  In  the  -  y  direction  with  the  dc  blae 
field  In  the  ♦  i  direction. 
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MISSION 
of 

Rome  Air  Development  Center 

RAVC  plank  and  zxecutet,  rzizan ch,  de.vzlopme.yvt,  tZkt  and 
selected  acquisition  programs  in  support  of  Command,  Control 
CcrmunecaUons  and  Intelligence  (C3 7}  activities.  Technical 
and  engineering  support  uiitfun  areas  of  technical  competence 
-c6  provided  to  ESD  Program  Offices  (POs )  and  other  ESV 
elements.  The  principal  technical,  mission  areas  are 
communication s,  electromagnetic  guidance  and  control,  Sur- 
veMancz  of  ground  and  aerospace  object s,  intelligence  data 
collection  and  handling,  information  system  technology, 
toted  state  sciences,  electromagnetics  and  electronic 
reliability,  maintainability  and  compatibility . 


